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ABSTRACT 

This  lecture  consists  of  three  parts:  1)  quantification  of  rates  of  gas-surface  interaction,  2)  formulation  of 
gas-surface  balance  conditions,  and  3)  survey  of  the  gas-surface  interaction  problems  in  the  entry  flights  of 
various  planets.  The  first  part  reviews  the  role  of  surface  rates  and  flow  parameters  in  Goulard 's  theory, 
the  catalytic  rates  for  Martian  entry  problem,  the  reaction  rates  for  oxidation  and  nitridation  of  carbon, 
rough,  surfaces,  and  the  relationship  between  forward  and  reverse  surface  rates.  The  second  part  covers  the 
method  of  expressing  the  diffusion  flux  at.  wall,  derivation  of  the  mass  balance  condition  at  wall,  and  its 
implementation  in  the  boundary  layer  and  computational-fluid-dynamics  formulations.  In  the  third  part.,  a 
brief  survey  is  made  of  the  gas-surface  interaction  problems  occurring  in  entry  flights  into  Earth,  Mars, 
Venus,  Titan,  and  outer  planets. 

1.0  NOMENCLATURE 

B  :  Rotational  constant,  cm-1, 

c  :  Speed  of  sound,  m/s. 

(c)  :  Crystalline  (solid)  state. 

G  :  Average  molecular  speed  (y/SkT/rrm),  m/s. 

Gp  :  Specific  heat  at  constant  pressure,  J/(kg-K). 

Cv  :  Specific  heat  at  constant  volume,  J /  (kg-K) . 

D  :  Effective  diffusion  coefficient. 

Da  :  Adsorption  energy,  J / mol. 

Df,  :  Diffusion  barrier  energy,  J/mol. 

Dg  :  Dissociation  energy  in  gas  phase,  J/mol. 

Dis  :  Multicomponent  diffusion  coefficient  between  species  i  and  s,  nr/sec. 

/  :  Stream  function,  Eq.  (32). 

Fs  :  A  bifurcation  parameter,  Eq.  (44),  dimensionless. 

g  :  Normalized  energy,  Eq.  (36). 

H  :  Enthalpy,  J/kg. 

AH  :  Energy  produced  by  reaction,  ev  or  J /mol. 

Js  :  Rate  of  mass  diffusion  of  species  s,  kg/(m2-sec). 

k  :  Boltzmann  constant,  J /K. 

kw  :  Catalytic  velocity,  Eq.  (5),  m/s. 

m  :  Mass  of  one  unspecified  particle,  kg. 

7 ns  :  Mass  of  one  particle  of  species  s,  kg. 

7iip  :  Pyrolysis  gas  mass  flow  rate,  kg/(m2-sec). 

7hptS  :  Mass  flow  rate  of  species  s  in  pyrolysis  gas,  kg/(m2-sec). 

7httS  :  Mass  flow  rate  of  species  s  resulting  from  reactions  and  sublimation,  kg/(m2-sec). 

7ht  :  Mass  flow  rate  of  material  removed  by  surface  reaction  and  sublimation,  kg/  (m2-sec) . 

7 hViS  :  Mass  flow  rate  of  species  s  in  sublimation  product,  kg/(m2-sec). 

Ms  :  Molecular  weight  of  species  s,  kg/mol. 

M  :  Average  molecular  weight,  kg/mol. 

Mw  :  Voo/cw,  dimensionless. 

7ia  :  Number  density  of  adsorbed  molecules,  in  2 . 

7 ig  :  Number  density  of  gas  molecules,  m~3. 

N  :  Number  of  species. 

p  :  Pressure,  Pascal. 

Pr  :  Prandtl  number,  dimensionless. 

q  :  Heat  transfer  rate,  W/m2. 
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Q  :  Partition  function. 

r  :  Radial  distance  from  axis  of  symmetry,  m. 

R  :  Nose  radius,  m. 

Re oo  :  Freestream  Reynolds  number  PooVooR/  poo- 

S(c)  :  Unspecified  solid  surface. 

Sc  :  Schmidt  number,  dimensionless. 

T  :  Temperature,  K. 

u  :  Tangential  velocity,  m/s. 

v  :  Normal  velocity,  m/s. 

Voo  :  Freestream  velocity,  m/s. 

x  :  Distance  along  wall  measured  from  the  stagnation  point,  m. 

Ws  ■  Rate  of  production  of  species  s,  sec  1 . 

Xs  :  Molar  fraction  of  species  s. 

y  :  Distance  normal  to  wall  measured  from  wall,  m. 

Zs  :  Virtual  species  mass  fraction  in  bifurcation  model,  Eq.  (47). 

as  :  Mass  fraction  of  species  s,  dimensionless. 

7  :  Gp/Cv. 

4>  :  Roughness  ratio  (BET  area/projected  area),  dimensionless. 

rj  :  Boundary  layer  coordinate  normal  to  wall,  Eq.  (31b),  dimensionless. 
r)a  :  Adsorption  (sticking)  coefficient,  dimensionless. 

rjc  :  Catalytic  recombination  coefficient,  dimensionless. 

//,.  :  Reaction  coefficient,  dimensionless. 

rfv  :  Evaporation  coefficient,  dimensionless. 

p.  :  Viscosity,  m2/sec. 

/ii  :  A  bifurcation  parameter,  Eq.  (48a),  dimensionless. 

p-2  ■  A  bifurcation  parameter,  Eq.  (48b),  dimensionless. 

p  :  Density,  kg/rn3. 

£  :  Boundary  layer  coordinate  along  wall,  Eq.  (31a),  dimensionless. 

1.1.  SUPERSCRIPT 

'  :  8/ dr). 

1.2.  SUBSCRIPTS 

oo  :  Freestream. 

a  :  Adsorbed  state, 

e  :  Edge  of  boundary  layer. 

E  :  Equilibrium. 

,  /  :  Forward  reaction. 

g  :  Gas. 

,  t  :  Reverse  reaction. 

R  :  Rotational  mode. 

T  :  Translational  mode. 

v  :  Evaporation. 

V  :  Vibrational  mode. 

w  :  Wall. 

2.0  PHYSICS  OF  GAS-SURFACE  INTERACTIONS 

2.1.  GOULARD’S  WORK  ON  CATALYTIC  EFFECTS 

The  problem  of  chemical  reactions  at  a  solid  surface  was  made  aware  first  by  Goulard.1  He 
pointed  out  that,  in  an  environment  where  chemical  reactions  are  frozen  in  the  boundary  layer,  heat 
transfer  rate  q  to  the  wall  can  be  reduced  by  making  it  noncatalytic  to  recombination.  In  most 
practical  flight  regimes  of  hypersonic  vehicles,  boundary  layer  flows  are  indeed  chemically  frozen. 
The  flight  experiments  conducted  with  the  Space  Shuttle  proved  Goulard’s  thesis. 

Goulard  considered  the  catalytic  recombination  process  typified  by  (  1  ev  =  96.48  J/rnol) 

0  +  0  +  S{c)  ~^02  +  5(e),  A H  =  5.08  ev  =  490  kJ/mol  (1) 
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in  the  stagnation  region  shown  schematically  in  Fig.  1.  Here,  S(c)  signifies  that  the  phenomenon 
occurs  at  the  surface  of  an  unspecified  material  (c  here  signifies  crystalline  state) .  His  work  introduced 
three  important  axioms.  The  first  axiom  is  that  the  ratio  of  heat  transfer  rates  between  a  fully 
catalytic  and  noneatalytic  surfaces  is  approximately  the  ratio  of  the  total,  i.e.,  translational  ( CpT ) 
plus  chemical,  enthalpy  to  the  translation-only  enthalpies. 

q  to  fully  catalytic  wall  CpT  +  energy  contained  in  dissociation  , 

q  to  noneatalytic  wall  CpT  v 


The  second  axiom  is  that  the  heat  transfer  rate  to  the  fully  catalytic  wall  is  approximately  the 
same  as  that  in  a  nonreacting  gas: 


q  to  a  fully  catalytic  wall  ps  q  in  a  nonreacting  gas. 


(3) 


The  third  axiom  is  that  the  transition  from  the  noneatalytic  to  catalytic  surface  occurs  at  around 
the  point  where  the  quantity 

 0.47  1 


Sc'2/3  pwku 

becomes  unity.  Here  k.w  is  the  so-called  catalytic  velocity  defined  as 


(4) 


kin  - 


VcJ 


a, 


rjcj  j  8kTw 
4  V  7T m 


0.34  r)Cifcw 


(5) 


where  r]cj  is  the  recombination  coefficient,  i.e.,  the  fraction  of  the  colliding  reactive  species  (e.g., 
atomic  oxygen)  that  undergoes  the  recombination  process.  The  subscript  “,/”  signifies  that  the 
coefficient  is  for  the  forward  process,  to  distinguish  it  from  the  reverse  process.  The  quantity  7  is 
the  specific  heat  ratio,  Cp/Cv,  and  c  is  the  sound  speed.  7  is  taken  to  be  1.4  to  arrive  at  the  factor 
0.34  in  Eq.  (5).  The  quantity  k.w  has  its  origin  in  the  kinetic  theory.  In  kinetic  theory,  the  number 
of  atoms  or  molecules  hitting  a  surface  is  given  by  the  so-called  arrival  rate 


.  n  n  8kT  _2 

Arrival  rate  =  —C  =  —\  -  m  sec 

4  4  V  7T771 


(6) 


where  n  is  the  number  density  in  in  3.  The  subscripts  e  and  w  in  Eq.  (4)  signify  the  edge  of  boundary 
layer  and  wall,  respectively. 

The  dependence  of  heat  transfer  rate  on  £  is 

Q  q,U,ncatalyUc  _  1 

Qfullycatalytic  diuincatalytic  1  C 

This  produces  a  curve  known  as  Goulard’s  S-curve,  which  is  shown  in  Fig.  2. 
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1000  100  10  1  0.1  0.01  0.001 

5  =  M/(r|Re0'5) 

Figure  2.  Goulard’s  S-curve  dependence  of  heat  transfer  rate  on 

By  assuming  that  the  pressure  distribution  in  the  stagnation  region  is  that  given  by  the  Newto¬ 
nian  hyersonic  theory,  the  quantity  due/dx  is  found  approximately  to  be 

dUe  Vqo  /r,PoQ  irj\ 

~dx  ~  R  V  1  ’ 

A  typical  value  of  f  can  be  found  by  taking  Sc  =  0.5,  Tc 0  =  250  K,  Te  =  6000  K,  Tw  =  3000  K, 
Poo/Pe  =  1/10,  p  oc  T3/i,  7  =  1.4,  and  defining  the  freestream  Reynolds  number 

_  Poo  boo  R 

U^oo  —  i 

P'00 


which  leads  to,  to  within  about  30%  accuracy, 


c 


1  Vpo  1 

Vcj  ('w  \/ R&oo 


(8) 


Writing  the  ratio  Voo/cw  as  Mw,  the  freestream  Mach  number  evaluated  at  the  wall  temperature, 
a  fully  catalytic  surface  and  a  noncatlaytic  wall  can  be  defined  as 

Fidly  catalytic  wall  :  r)cj  %/Rb0 0  >>  Mw  (8 a) 

Noncutalytic  wall  :  r]cj  \/Re oo  <<  Mw.  (86) 

Using  arc-heated  wind  tunnels,  rjcj  for  the  0  +  0—^  Oa  and  N  +  N  — >■  N2  reactions  have  been 
measured  for  several  common  heatshield  materials.  Generally,  it  is  between  0.1  and  1  for  a  metal, 
between  0.01  and  0.1  for  a  metal  oxide,  and  less  than  0.01  for  silica  and  alumina,  two  most  commonly 
found  glassy  oxides. 

It  is  to  be  noted  here  that  the  surface  of  most  innoble  metals,  such  as  copper,  turns  quickly 
into  a  metal  oxide  when  it  is  exposed  to  a  hot  stream  containing  atomic  oxygen.  A  systematic 
study  was  made  in  Ref.  2  to  find  out  how  fast  a  copper  surface  turns  into  a  copper  oxide  surface. 
The  experiment  was  made  at  a  pressure  of  the  order  of  0.01  atm.  Within  1  sec  at  this  pressure, 
copper  became  copper  oxide  when  it  was  tested  in  an  arc-heated  wind  tunnel.  In  Ref.  3,  copper 
calorimeter  was  used  to  measure  the  heat  transfer  rate  at  a  pressure  of  2/3  atmospheres.  In  this 
case,  enthalpy  of  the  flow  was  determined  spectroscopically.  The  measured  heat  transfer  rate  to  the 
copper  calorimeter  was  consistent  with  the  assumption  that  .  is  0.01,  which  is  a  typical  value  for 
metal  oxides. 

Of  particular  interest  is  the  r]cj  value  for  O-atoms  on  the  surface  of  graphitic  carbon,  which  is 
used  as  an  ablative  heatshield.  Efforts  have  been  made  in  the  past  to  measure  this  quantity  in  an 
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atomic  beam  experiment  (see,  e.g.,  Ref.  4).  Those  experiments  found  no  O2  in  the  product:  the 
O-atoms  either  formed  CO  by  combining  with  carbon  or  reflected  elastically  as  an  O-atom,  without 
forming  O2.  That  is, 

For  O  +  O  +  C(c)  — >  0-2  +  C(c)  :  r]cj^0.  (9) 

2.2.  r)cj  FOR  CO  +  O  — >  C02 

When  unmanned  Mars  entry  missions  were  started,  question  was  raised  as  to  what  is  the  impact 
of  surface  catalytic  recombination.  Martian  atmosphere  consists  mostly  of  CO2.  At  the  edge  of 
boundary  layer,  CO2  is  fully  dissociated  into  CO  and  O.  CO  and  O  could  conceivably  recombine  as 

CO  +  O  +  5(c)  — >  CO-2  +  5(c)  +  5.35  ev. 

At  a  flight  velocity  of  6  km/s,  a  typical  entry  velocity  of  a  Martian  vehicle,  the  kinetic  energy  of 
one  CO2  molecule  is  8.21  ev.  Therefore  the  reaction  energy  5.35  ev  is  more  than  half  of  the  flow 
enthalpy.  If  CO  and  O  do  not  recombine  catalytically  at  wall,  Eq.  (1)  tells  us  that  the  heat  transfer 
rate  will  be  less  than  half  of  the  fully  catalytic  case.  That  is,  reward  to  making  the  wall  noncatalytic 
is  quite  large.  The  same  problem  exists  for  Venusian  entry  because  the  atmospheric  composition  of 
Venus  is  similar  to  that  of  Mars. 

Automobile  manufacturers  have  developed  catalytic  converters  which  burn  CO  with  O2  to  make 
C02-  They  accomplished  it  by  making  CO  molecules  collide  many  times  with  platinum  surfaces  in 
the  presence  of  02-  This  means  that  the  catalytic  efficiency  r]cj  is  very  small  for  CO  +  O2.  NASA 
Ames  Research  Center  carried  out  an  experiment  to  measure  rjcj  for  the  CO  +  O  reaction.  Because 
the  surface  of  the  heat-shield  for  Martian  entry  vehicle  is  usually  fused  silica,  the  experiment  was 
conducted  for  the  silica  surface.'"'  The  experiment  found  that  the  reaction  CO  +  O  is  too  slow  to  be 
detected. 

That  r)cj  for  CO  +  O  — >■  CO2  is  small  can  be  deduced  theoretically.  A  catalytic  reaction  involves: 
1)  adsorption  (sticking)  of  the  two  component  species,  2)  migration  of  the  component  species  in  the 
adsorbed  phase,  3)  recombination  of  the  two  partners,  and  4)  desorption  of  the  product.  The  CO 
and  O  particles  that  strike  the  wall  may  not  be  adsorbed,  i.e.,  fail  step  1),  or,  skip  step  3)  and  desorb 
without  recombining.  The  rates  of  these  four  processes  can  be  modeled.  That  work  will  be  done  by 
author  lecturers.  Here,  only  the  steps  1)  and  4)  will  be  given  to  explain  why  r]cj  for  the  process  is 
slow. 

The  r]cj  for  the  catalytic  process  cannot  be  larger  than  the  adsorption  coefficient  of  CO,  because 
the  latter  is  a  component  of  the  former.  For  a  monatomic  gas,  the  adsorption  coefficient,  denoted 
here  rja ,  is  known  to  be  a  function  of  the  ratio  of  the  mass  of  the  gas  atom  mg  to  that  of  the  atoms 
forming  the  solid  mc.  Classical  calculation  of  the  surface  processes6  predict  that,  when  the  mass 
ratio  is  much  smaller  than  unity,  r)a  behaves  as 


The  existing  experimental  data  with  helium,  neon,  and  argon7  agree  roughly  with  this  theory. 

As  mentioned,  the  surface  of  the  heatshield  for  the  Martian  entry  vehicles  is  usually  Si02- 
Because  the  molecular  weight  of  Si  is  28,  collisions  of  the  CO  and  O  will  have  little  effect  of  this 
mass  ratio.  Thus,  in  the  absence  of  other  effects,  i]a  should  be  nearly  1  in  the  Martian  case. 

But  there  is  what  is  known  as  entropy  effect,  which  concerns  the  change  in  the  internal  degrees 
of  freedom  during  adsorption.  (To  change  internal  degrees  of  freedom,  heat  must  be  given  or  taken 
away  at  a  constant  temperature.  By  the  definition  of  entropy,  this  implies  a  change  in  entropy.) 
Let  us  designate  the  gas  phase  quantities  by  a  subscript  g  and  the  adsorbed  state  quantities  by  a 
subscript  a.  The  adsorbed  gas  species  are  held  there  by  the  attractive  potential  between  the  solid 
surface  and  the  gas,  Da.  The  numerical  value  of  Da  is  known  for  some  cases8.  For  CO  on  a  metal, 
it  is  typically  between  1.5  and  6  ev. 

A  molecule  in  the  adsorbed  state  oscillates  within  the  adsorption  potential  well  approximately 
with  the  same  frequency  as  the  solid  lattice,  which  is  v  k,  k,Tw/h  on  the  average.  The  fraction  of 
the  adsorbed  CO  molecules  that  has  an  energy  exceeding  the  adsorption  energy  Da  is  exp(-£>a/A:Tl„). 
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Therefore,  the  number  of  CO  molecules  in  the  adsorbed  state  that  desorb  from  the  adsorption 
potential  well  within  1  sec  per  1  m2  is 

Desorption  rate  =  ™  exp(—  °  )na  m~2sec~1 

lb  f\j-L  uj 

As  given  by  Eq.  (6),  the  flux  of  gas- phase  molecules  reaching  the  wall  is  Gwiig/ 4.  Only  the  r]a 
fraction  of  this  flux  will  be  adsorbed,  so  that 


Adsorption  rate  =  -j- Cwng  in  2 sec  1. 


Under  equilibrium,  these  two  rates  must  be  equal,  i.e., 


which  leads  to 


i  I  'w  /  Da  !ja  ( 

— 1 =  — ’*» 


4  h  Z  \r  Da  na 

da  =  — -^ea:p(-— )  — . 


Gw  h 


id  \r  ng 


The  ratio  n^g/flo  is  given  by  the  partition  functions  for  translation  Qt,  for  vibration  Qv>  and 
for  rotation  Qr,  by 

naE  QTaQvaQlla  / 

-exP(u ir)- 


Therefore 


ng  ()  j  (I  ()  \  ' g  ()  jin  i  i  'u 
4  kTw  i)  j  V,  C)  \  'a  () lia 

G-w  h  ()■]  g  ()  Qun 


da  = 


Here,  all  partition  functions  are  evaluated  at  the  wall  temperature  Tw. 

The  translational  partition  function  in  the  gas  phase,  Qxg,  is  given  by 

QTg={2-^kf/  2 

The  largest  possible  value  of  the  translational  partition  function  in  the  adsorbed  phase,  Qtu,  occurs 
when  the  adsorbed  molecule  has  a  full  freedom  of  translation  in  two  dimensions,  in  which  case 


Qraifull  2D)  = 


2  mnkTw 
h? 


When  there  is  an  energy  barrier  of  Df,  that  hinders  the  two-dimensional  motion,  so  that  only  those 
molecules  with  energy  greater  than  Df,  have  the  two-dimensional  translational  freedom,  one  must 
multiply  a  factor  exp(-Df,/kT)  to  this  value,  and  therefore 


27T  mkTw  Db 

Q^  =  ~ M?- 


(10) 


This  leads  to 


4  kTw  ,  h  \l/2  QvaQ Ra  /  \ 

Cw  h  ^27rmkTw'  Qv9Qr9 ^ 


QvaQ  Ra  /  4^5  , 

QvgQ^-^J 


(11) 


For  a  monatomic  gas,  Qv  and  Qr  are  both  unity.  If  Db  =  0,  Eq.  (11)  gives  r)a  =  1.  This  agrees 
with  what  is  known  about  r)a  for  monatomic  species. 

At  best,  the  rotational  motion  of  an  adsorbed  molecule  occurs  around  only  one  axis,  the  axis 
normal  to  the  surface.  By  assuming  that  the  vibrational  motion  of  the  adsorbed  phase  is  excited  to 
the  same  extent  as  the  gas  phase,  one  obtains  the  largest  possible  value  of  the  adsorption  coefficient 
to  be  _ 


da 


Db_ 

kT,„ 


(12) 
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Assuming  Da  =  0,  the  numerical  value  of  pa  for  CO  at  the  wall  temperature  of  2000  K  is  deduced 
to  be  0.0373.  If  Da  is  finite,  or  the  rotational  motion  is  hindered  completely,  or  if  the  vibrational 
motion  is  hindered  in  the  adsorbed  phase,  r)a  will  be  smaller.  This  explains  at  least  partly  why  the 
catalytic  recombination  coefficient  for  CO  +  O  is  small  on  the  surface  of  Si02- 


2.3.  GAS-SOLID  REACTIONS  AND  SUBLIMATION 


At  high  heating  rates,  the  heatshield  surfaces  ablate.  For  an  ablating  environment,  carbonaceous 
heatshield  is  commonly  used.  For  this  case,  the  gas  species  O,  N,  and  H  can  combine  with  the  carbon 
atoms  of  the  solid  surface  as 


O  +  C(c)  -»•  CO, 

AH  =  3.74  ev, 

(13) 

N  +  (7(e)  -»•  CN, 

AH  =  0.35  ev, 

(14) 

H  +  (7(c)  -»•  CH, 

AH  =  -3.88  ev, 

(15) 

H  +  2 (7(c)  -»•  C2H, 

AH  =  —2.67  ev. 

(16) 

For  reaction  (13),  the  existing  experimental  data9  show  that  the  reaction  coefficient  7)rj  can  be 
approximated  by 


O  +  (7(c)  -►  CO  :  rjrJ  «  0.63erp(-1160/T»,). 
For  reaction  (14),  Ref.  10  shows 


N  +  C(c)  -»•  CN  :  7i,. j  «  0.3. 

Reactions  (15)  and  (16)  can  occur  only  at  high  wall  temperatures,  and  have  not  yet  been  studied. 
However,  the  maximum  possible  values  can  be  reasoned  simply  from  the  requirement  that  only  the 
particles  with  the  required  energies  (3.88  ev  for  H  +  C(c)  and  2.67  ev  for  H  +  2C(c))  can  produce 
the  change.  This  leads  to 


H  +  (7(c)  -»•  CH  :  7]rJ  <  exp ( — 45 , 000 /Tw ) 

H  +  2 (7(c)  -»•  C2H  :  r)rJ  <  exp(-?A,  000 /Tw). 

If  the  flow  contains  CH  or  C2H  at  the  boundary  layer  edge,  as  in  the  atmosphere  of  Titan,  reactions 
(15)  and  (16)  will  occur  in  the  reverse  direction  at  wall.  That  is,  carbon  will  condense  on  the  wall 
from  CH  and  C2H. 

At  high  wall  temperatures,  evaporation  (sublimation  if  the  surface  is  solid)  occurs.  Imagine 
that  a  wall  material  is  enclosed  in  a  high  temperature  furnace  and  is  heated.  Evaporation  will 
occur  and  fill  the  inside  volume  with  the  vapor  of  the  material.  The  equilibrium  vapor  pressure  for 
most  materials  are  known,  and  are  given  in  the  JANAF  Thermochemical  Tables.11  Denoting  the 
equilibrium  vapor  number  density  by  n_g,  the  flux  of  the  vapor  species  hitting  the  wall  is  given  by 
the  kinetic  theory  to  be  HeCw/ 4.  Of  this  flux,  in  general  an  rfv  fraction  will  be  adsorbed  on  the  wall. 
Under  equilibrium,  this  flux  must  equal  the  flux  of  evaporation.  Thus  the  evaporation  rate  is  given 
by  b  Under  nonequilibrium,  the  vapor  number  density  n  will  produce  adsorption  rate  of 

rjrn(  'w  1.  Therefore,  the  net  mass  flux  of  evaporation  will  be 

Evaporating  flux  =  7-^-Gw(nE  —  n)  m_2aec_1 

where  ng  and  n  are  evaluated  at  the  wall.  The  corresponding  mass  flux  becomes,  in  terms  of  mass 
fraction  a, 

J  =  -^PwCw(aE  -  a)  kg/(m 2  -  sec )  (17) 

Graphitic  carbon  vaporizes  into  C,  C2,  O3,  C4,  and  C5.  Of  these,  the  predominant  reaction  is 

3 (7(c)  -»•  (7a,  A H  =  -7.51  ev  (18) 
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For  this  process,  i]r  has  been  measured  in  Ref.  12.  The  results  can  be  fitted  by 

Vv(C3)  =  30ea;p(-21490 /Tw). 

Glassy  oxides  liquefy  before  they  vaporize.  Vaporization  of  liquid  oxides  involves  at  least  two 
steps.  For  instance,  when  Si02  is  heated,  its  vapor  consists  of  SiO  and  O2.  No  quantitative  infor¬ 
mation  is  available  for  r/v  for  this  class  of  materials. 

2.4.  ROUGH  SURFCES 

When  the  surface  reaction  of  the  type  Eqs.  (13)  to  (16)  or  vaporization,  Eq.  (18),  occurs, 
the  surface  becomes  rough.  The  extent  of  roughness  is  measured  by  the  ratio  of  the  so-called  BET 
(Bennett-Emmons- Thomas)  area  to  the  projected  area,  named  here  4> 

BET  area 

Roughness  ratio  cp  =  - - - .  (19) 

projected  area 

BET  area  is  the  wet  area  of  a  substance,  and  is  determined  by  an  experiment  shown  schematically 
in  Fig.  3.  One  places  the  sample  inside  a  sealed  jar  and  evacuates  the  jar.  The  sample  is  heated 
while  pumping  to  purge  all  gas  molecules  adsorbed  on  it.  Then  one  admits  gas,  typically  nitrogen, 
while  measuring  the  mass  flow  rate  that  enters  the  jar  and  watching  the  pressure  inside  the  jar.  In 
the  beginning,  pressure  in  the  jar  does  not  rise  because  all  gas  molecules  are  consumed  in  adsorbing 
on  the  wet  surface  of  the  sample.  After  a  while,  pressure  begins  to  rise,  because  the  surface  is  fully 
covered  by  a  monolayer  of  the  adsorbed  gas.  The  surface  area  occupied  by  one  molecule  is  typically 
10-15  cm2.  By  dividing  the  total  number  of  gas  molecules  consumed  up  to  this  point  by  this  value, 
one  obtains  the  BET  area. 


GAS 


\MOUM  OI;  GAS  ADMIT!  TD  INTO  JAR 


Figure  3.  Schematic  of  the  procedure  to  determine  the  BET  area  of  a  solid  surface. 

The  BET  area  has  never  been  measured  during  ablation.  After  ablation,  theoretically  it  can  be 
measured.  However,  a  heatshield  material  becomes  highly  porous  after  ablation.  The  experiment 
shown  in  Fig.  3  will  measure  the  wet  area  of  the  pores  as  well  as  the  area  exposed  to  the  boundary 
layer  flow.  Therefore,  the  roughness  ratio  4>  can  only  be  guessed.  For  many  such  porous  materials, 
4>  is  of  the  order  of  100.  Perhaps  we  should  consider  <f>  to  be  between  10  and  100. 

We  can  now  conduct  a  thought  experiment  as  to  what  will  be  the  magnitude  of  the  adsorption 
coefficient  r)a  on  a  rough  surface,  as  shown  schematically  in  Fig.  4.  Here,  a  swarm  of  atoms  is  hitting 
the  inside  of  a  crevice.  In  the  first  collision,  r)a  fraction  is  adsorbed,  and  1  -  rja  fraction  leaves.  This 
remaining  swarm  of  atoms  hits  the  wall  the  second  time.  The  atoms  that  leave  are  now  (1  -  r)a)2 
fraction  of  the  original  swarm.  After  n  number  of  collisions,  there  will  be  (1  -  r)a)n  fraction  of  the 
original  swarm  left  unadsorbed.  Therefore,  the  fraction  of  the  atoms  adsorbed  are 

Effective  adsoi'ption  coef  =  1  —  (1  —  rja)n.  (20) 

The  total  number  of  collisions,  n,  will  be  of  the  order  of  <f>.  So,  to  this  order  of  accuracy,  n  can  be 
replaced  by  0.  If  r)a  is  very  small  and  if  f>  is  large,  Eq.  (20)  can  be  approximated  by 

Effective  adsorption  coef  «1-  exp{—r)a4>) .  (21) 
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These  show  that,  when  a  surface  is  rough,  its  effective  r)a  approaches  unity. 


Figure  4.  Schematic  showing  collisions  of  particles  on  a  rough  surface. 

2.5.  INTER-RELATIONSHIP  AMONG  r/’s 

The  above  reasoning  can  be  applied  to  catalytic  reactions  and  surface  reactions  to  conclude  that 
not  only  r)a,  but  rjc ,  i]r .  and  rj,:  all  approach  unity  when  roughness  ratio  is  high.  However,  all  ?fs 
cannot  become  unity  together  because,  by  definition,  the  sum  of  all  rf s  for  one  species  must  not  be 
larger  than  unity: 

£77  <  1.  (22) 

Otherwise,  mass  conservation  is  violated. 

Take  an  example  with  O-atoms  hitting  a  carbon  surface.  O-atoms  can  possibly  recombine  on  the 
surface  with  a  coefficient  r]cj  and  simultaneously  can  react  with  carbon  to  form  CO  with  a  coefficient 
T\rJ.  The  sum,  Tjej  +  rjrj,  cannot  be  larger  than  unity.  Because  rjcj  is  nearly  zero  and  r)rj  is  large, 
the  best  practice  may  be  to  take  the  effective  ijrj  to  be  unity  and  the  effective  r]cj  to  be  0  for  a 
rough  carbon  surface. 


2.6.  FORWARD  VS.  REVERSE  RATES 

In  Eq.  (17),  the  two  terms  in  the  parenthesis  represent  the  forward  and  reverse  rates  in  evapo¬ 
ration.  Likewise,  for  both  the  catalytic  recombination  processes  typified  by  Eq.  (1)  and  the  surface 
reactions,  Eqs.  (13)  to  (16),  the  reverse  rates  must  be  accounted  for.  The  reverse  rates  are  chosen  so 
as  to  satisfy  the  principle  of  detailed  balance  in  equilibrium.  Otherwise,  equilibrium  condition  will 
not  be  reached  when  it  should  be.  This  aspect  of  the  problem  is  discussed  in  Ref.  13. 

For  example,  let  us  take  the  recombination  process,  Eq.  (1).  There  is  an  opposing  (reverse) 
reaction 

0  +  0  +  S(c)  02  +  S(c).  (23) 

The  fraction  of  the  O2  molecules  that  undergoes  this  process  will  be  designated  by  //Ct,. .  The  net 
reaction  rate  will  be  the  difference  between  the  forward  and  the  reverse  rates 


Net  rate  of  0  removal  =  ^^-CWtono  ~  2~^-&v},o3nOt 


Here,  the  factor  2  in  the  second  term  signifies  that  two  O-atoms  are  created  by  the  dissociation  of 
one  O2  molecule.  The  mass  flux  of  the  O-atms  removed  by  this  process  can  be  written  as 

Jo  —  hy  >OOtw  w) .  (24) 


Under  equilibrium,  the  forward  and  reverse  rates  must  be  equal,  which  leads  to 

Vcj  Cw  o  no,E  1  ,  no  . 

Vo,r  =  - " -  =  -j=Tktf(  —  )E. 

O 2  ^02 , E  A/"  ^02 


(25) 
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The  ratio  no  EftlOsM  is  given  by  the  equilibrium  relationship 


( Qt,o )2 


Qt,o-2Qv,o2Q  r,o2 


e®P(- 777H 

m  7/i 


from  which  there  results 

(72^)e  =  1  (Qt,o)2  exp(  Dg  ) 

no-2  no  Qi  a  >2Q\\o-,  Qu,o-2  kTw 


A  similar  relation  holds  for  all  other  catalytic  recombination  processes. 
Similarly,  for  the  surface  reaction  (13),  the  net  rate  is 


Net  rate  of  O  removal  =  ^-^-Cw,ono  —  CWtconco  in  2sec  1 


Here 


Jo  —  ^Pw(jlr,f(Jw,0(X0,w  Vr,r^J'W,CO^GO,w)  k‘9  /  (in  Sec). 

(J-w.o  /  no  \ 

Vr,r  —  VrJ-p^  (~  )e- 


CWtco  nco 

The  ratio  (no /nco) e  is  given  by  the  equilibrium  relationship 


=  1  Qt,oQt,c  Dg  ^ 

nco  nc,E  Qt,coQv,coQr,co  MTw 


(26) 


(27) 

(28) 


(29) 


where  jic,e  is  the  equilibrium  vapor  number  density  of  Cg.  A  similar  relationship  holds  for  all  other 
surface  reactions. 


2.7.  HEAT  TRANSFER  RATES 

All  three  surface  phenomena,  i.e.  catalytic  recombination,  reaction,  and  sublimation,  produce 
heat  corresponding  to  AH.  All  these  heat  enter  into  the  heatshield  wall.  The  heat  transfer  rate 
caused  by  these  phenomena  are 


q  =  Rate  of  reaction  in 


event  AH  ■m  J/mol 
m2sec  X  6.0225  x  1023 


W/m 2. 


(30) 


Here,  one  urns  not  confuse  event/(m2sec)  with  rate  of  species  removal  or  species  production.  For 
instance,  for  the  oxygen  recombination,  Eq.  (9),  removal  of  two  atomic  oxygen  equals  one  event. 
Therefore,  even!  (m2sec)  is  half  of  the  rate  of  removal  of  O-atoms.  The  above  equation  is  valid 
unconditionally,  regardless  of  the  magnitude  of  ablation  rate. 


2.8.  EXTENTION  OF  GOULARD’S  THIRD  AXIOM 

We  can  now  ponder  how  much  impact  the  coefficients  1],.  (in  the  surface  reactions  of  the  form 
(13)  to  (16))  and  //,  (in  sublimation  of  the  form  (18))  will  have  on  heat  transfer  at  the  stagnation 
point.  In  other  word,  we  would  like  to  extend  Goulard’s  analysis  to  surface  reactions  and  sublima¬ 
tion.  Unlike  catalytic  recombination  process,  surface  reaction  and  sublimation  produce  ablation.  To 
extend  Goulard’s  theory  to  an  ablating  surface,  we  have  to  start  with  the  fundamental  conservation 
equations. 

Fundamental  conservation  equations  for  the  frozen  boundary  layer  flow  are  given  in  Goulard’s 
work1,  and  so  will  not  be  repeated  here.  Those  equations  are  valid  even  when  surface  reactions  and 
sublimation  occur.  Following  Goulard,  we  write  those  conservation  equations  using  the  coordinates 
x  and  y  shown  in  Fig.  1.  Then  one  introduces  the  coordinate  transformation 


£.=  I  Petieiier2dx, 

Jo 


(31a,  b) 
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and  the  stream  function  /  that  satisfies  the  conditions 


3$ 


—  =  pur,  —  =  - pvr 

dy  ox 


The  stream  function  /  can  be  written  as 

/  =  - 


pv 


y/2pepe(due/dx) 


(32) 


(33) 


Then,  denoting  the  partial  differentiation  d/d rj  by  ',  the  mass  and  momentum  conservation  equations 
reduce  to: 


(— r)'+/r+k--/,2)  =  o. 

PePe  2  p 

Following  Goulard,  we  assume  pe/p  =  f'2,  so  that  Eq.  (34)  reduces  further  to 


+  //"  =  0. 

PePe 


(34) 


(35) 


The  energy  equation  is  transformed  likewise  by  introducing  a  normalized  energy  function 

g  =  H/He  (36) 


into 


{—^g'Y + fa* d4l(i-  =  o. 

PepePr  He  Pr  pep.e 


Assuming  Prandtl  number  to  be  unity,  Eq.  (37)  reduces  further  to 

1  a'Y  +  fa'  =  o. 


PePe  Pr 

The  species  equation  is  transformed  likewise  into 


PP  1  \  /  i  j*  /  n 

:«s)  +  Jas  —  6 


y  PePe  Sc 


(37) 


(38) 


(39) 


When  there  is  no  ablation,  the  wall  value  of  /is  zero  according  to  Eq.  (33) .  The  numerical  value 
of  0.47  appearing  in  Goulard’s  parameter  f,  Eq.  (4),  is  the  wall  slope  of  species  mass  fraction  a  in 
the  solution  of  Eq  (39),  ( dot/ dr/) w ,  for  Sc  =  1.  0.47  also  happens  to  be  the  wall  slope  of  the  energy 
variable  in  Eq.  (38),  ( dg/drj)w ,  for  Pr  =  1. 

When  there  is  ablation,  the  wall  value  of  / 


fw 


Pw'Vw 

^/2pepe(due/dx) 


(40) 


is  no  longer  zero,  but  is  finite  negative.  When  fw  has  a  finite  negative  value,  numerical  value  of 
{da/drf).w  and  (dg/drj).w  are  different  from  0.47.  Numerical  solution  of  those  differential  equations 
give  these  wall  slope  values  as  shown  in  Fig.  5.  A  two-parameter  lit 


f9g_,  _ 

{dy)w'  {dph 


0A7ewp(Q.mfw  -  4.017 /2 ) 


(41) 


is  shown  in  the  figure  for  comparison.  The  analytical  fit  is  seen  to  be  valid  to  fw  of  up  to  about  -1.2. 
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Figure  5.  Wall  slope  of  energy  and  species  equations  at  the  stagnation  point, 

under  Goulard’s  simplifications. 

The  wall  slope  of  the  energy  equation,  (dg/drj)w,  is  a  measure  of  heat  transfer  rate.  As  Fig. 
5  shows,  it  decreases  rapidly  as  -fw  increases.  This  phenomenon  is  known  as  convective  blockage 
effect.  Using  Eq.(7),  fw  can  be  written  as 


U  = - 

Poo  v  co  Pe 


where  m  is  the  rate  of  mass  injection  by  ablation.  Taking  Poo/pe  to  be  0.1,  this  becomes 


777  - 

/„»-0.33  — —  y/feZ.  (42) 

Poo  *  CO 

The  quantity  fw  is  the  fundamental  quantity  in  the  so-called  blowing  parameter,  which  is  a  ratio  of 
the  mass  injection  rate  m  and  the  freestream  flow  rate  PooVoo  and  which  is  defined  slightly  differently 
by  various  researchers. 

By  replacing  0.47  in  Goulard’s  £,  Eq.  (4),  by  the  expression  (41),  one  has  Goulard’s  third  axiom 
extended  to  the  case  of  finite  ablation  rate.  Generally,  finite  fw  pushes  the  S-curve  in  Fig.  2  toward 
the  right,  so  that  the  wall  behaves  more  as  a  catalytic  wall  for  the  same  catalytic  coefficient. 

The  S-curve  can  be  used  also  to  gage  the  impacts  of  surface  reactions  or  sublimation  on  heat 
transfer.  For  instance,  if  there  was  surface  reaction  of  the  form  Eq.  (13)  but  no  surface  catalysis  or 
sublimation,  heat  transfer  rate  will  increase  due  to  the  surface  reaction.  The  S-curve  in  Fig.  2  is 
applicable  exactly  as  it  is.  However,  the  magnitude  of  the  heat  transfer  rate  will  have  to  be  changed 
using  the  curve  shown  in  Fig.  5. 


3.0.  KINETIC  BOUNDARY  CONDITIONS 
3.1.  BIFURCATION  MODEL 

When  the  heating  rate  is  above  about  500  kW/rn2,  a  heff shield  material  loses  mass  during 
heating.  That  is,  it  is  not  truly  reusable.  Between  about  500  kW/rn2  and  1.5  MW/rn2,  the  mass  loss 
rate  is  so  small  that  the  material  can  be  replenished  after  use,  to  be  used  the  next  time.  When  the 
heating  rate  exceeds  about  1.5  MW /nr2,  the  mass  loss  rate  is  so  high  that  an  ablating  material  must 
be  used.  An  ablating  material  consists  generally  of  two  components,  a  highly  temperature-resistant, 
highly  tensile-resistant  fibre  matrix  and  a  resin.  The  material  is  usually  porous.  When  it  is  heated, 
resin  vaporizes  at  relatively  low  temperatures,  e.g.,  500  K.  The  vapor  passes  through  the  pores  and 
escapes  from  the  surface  (see  Fig.  6).  The  vaporization  process  absorbs  heat.  While  passing  through 
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the  pores,  the  vapor  molecule  decomposes  into  smaller  molecules.  This  decomposition  process  absorbs 
heat  and  thereby  cools  the  material  farther.  This  gas  is  called  pyrolysis  gas,  and  the  process  is  called 
pyrolysis  process.  In  the  region  near  the  wall  surface,  the  matrix  is  devoid  of  the  resin.  This  region 
is  called  char.  The  zone  where  the  vaporization  of  the  resin  is  occurring  is  called  the  pyrolysis  zone. 


Diffusion  of  O 


Figure  6.  Schematic  of  the  phenomena  occurring  below,  at,  and  above  the  wall  surface  of 

an  ablating  heatshield. 

The  general  principle  of  mass  balance  at  wall  states  that  whatever  is  produced  or  removed  at 
wall  is  taken  away  or  supplied  by  diffusion  in  the  boundary  layer,  i.e. 


Rate  of  production  at  wall  =  rate  of  removal  by  dif fusion. 


(43) 


Therefore,  before  proceeding  any  further,  we  must  know  how  to  express  the  rate  of  diffusion. 

The  rate  of  diffusion  of  a  species  is  given  in  general  by  the  Chapman-Enskog  theory  using 
the  multi-component  diffusion  coefficients.  According  to  this  formulation,  the  rate  of  diffusion  of 
species  s  is  a  weighted  sum  of  gradients,  dXi/dy,  of  all  species  except  s.  The  weights  contain  the 
multi-component  diffusion  coefficients  D/s  as  a  factor.  This  complicates  wall  mass  balance  condition 
greatly. 

A  considerable  simplification  is  achieved  by  adopting  the  so-called  bifurcation  model.  The  model 
is  based  on  the  assumption  that  the  multipcomponent  diffusion  coefficient  between  species  i  and  j 
can  be  approximated  by 


Dt 


D 

TTf 


(44) 


Here,  D  is  a  reference  diffusion  coefficient,  which  is  chosen  in  Ref.  14  to  be  the  self  diffusion  coefficient 
of  N2.  Eq.  (44)  is  theoretically  rigorous  if  the  interaction  among  species  is  purely  by  the  induced 
dipole- induced  dipole  interactions.  At  temperatures  of  up  to  about  3000  K  in  air,  molecular  potentials 
are  indeed  almost  those  of  induced  dipole-induced  dipole  interactions.  At  higher  temperatures,  that 
is  not  so,  and  so  a  very  rough  approximation  must  be  made  with  Di  /s  in  order  to  satisfy  this 
condition.  The  values  of  Ffs  so  selected  are  given  in  Ref.  14.  Fs  is  proportional  approximately  to 
the  square- root  of  its  molecular  weight  relative  to  the  reference  species  (which  is  0.028  kg/mol  for 

N2): 


Fk 


Ms 

0.028' 


(45) 


If  Eq.  (44)  holds,  the  bifurcation  model  states  that  the  rate  of  mass  diffusion  of  species  s  can 
be  written  as 

pFtp'2  d%s  t//2\  { iC\ 

kg /(in  (46) 


J,= 
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where  Z.  is  a  virtual  mass  fraction 


and 


M  is  the  average  molecular  weight 


^s  = 

M8xs 

Fsp i  ’ 

(47) 

Mi  —  ^  ^  XSF 5, 

v—'  MSXS 

F2  =  2^  f  • 

1  s 

(48a,  b) 

M  =  J2xbM8. 


The  rate  of  production  of  species  s  by  chemical  reaction  will  be  denoted  by  W8.  The  conservation 
equation  for  species  s  can  be  written  as 

das  das  d  ,  pDp ,3  dZs 


^  dx  P1  dy  dyx  Mfj-i  dy 


■)  +  pWs. 


The  quantity  Dp%fiMp\)  can  be  considered  to  be  an  effective  diffusion  coefficient: 

Df.l-2 


D  = 


Mp  i 


(49) 

(50) 


Schmidt  number  can  now  be  defined  as 

Sc  =  p/pD.  (51) 

Sc  is  a  function  of  temperature,  but  is  not  a  function  of  species.  The  virtual  mass  fraction  Zs  in  Eq. 
(47)  can  be  written  further  as 

ry  _  M8X8  Yhk  MkXk  _  X k  MkXk  _  M 

S  ~  TkMkXk  Fsp-2  ~ as  Fsp-2  ~ a 
One  can  introduce  a  parameter  C8: 


Fsp-2 


Cs  = 


M 

FsP-2 


Using  these  notations,  the  species  conservation  equation,  Eq.  (49),  can  be  written  as 

da3  dots  d  r  p  d{Csots ), 

pu-F  +  po  — £  =  -[±-  s  +  pWs. 

dx  dy  dy  Sc  dy 


(52) 

(53) 


3.2.  MASS  BALANCE 

In  order  to  describe  the  behavior  of  an  ablating  wall,  what  happens  inside  the  heatsliield  material, 
i.e.,  material  response,  needs  to  be  calculated.  This  can  be  done  by  using  a  computer  code  such  as 
Super  Charring  Materials  Ablation  (SCMA)  code  given  in  Ref.  15.  The  calculation  will  yield  the 
mass  flow  rate  of  the  pyrolysis  gas  issuing  from  the  wall,  mp.  By  invoking  equilibrium  among  the 
various  species  constituting  the  pyrolysis  gas,  one  can  determine  the  mass  flow  rates  of  each  of  the 
species  mPjS.  The  mass  fraction  of  species  s  in  the  pyrolysis  gas  will  be  represented  by  a8tP. 

The  total  mass  flow  rate  emerging  from  the  wall  is  the  sum  of  the  pyrolysis  mass  flow  rate  mv 
and  the  mass  flow  rate  produced  at  wall  m* 

PwVw  =  mp  +  jht,  (54) 

The  surface  mass  removal  rate  mt  is 

mt  =  sum  of  mass  removal  rates  by  surface  reaction 
+sum  of  mass  removal  rates  by  sublimation 

As  can  be  seen  from  Eqs.  (24)  and  (27),  mass  removal  rates  by  surface  reactions  consist  of  terms 
proportional  to  a8  .w.  As  can  be  seen  from  Eq.  (17),  the  mass  removal  rates  by  sublimation  consist 
of  terms  proportional  to  aSjW  and  terms  independent  of  aSiW.  Thus,  one  can  write 

m  =  y~$58as,w  +  es).  (55) 

S 

The  parameters  Ss  and  es  are  determined  by  adding  up  contributions  from  all  surface  chemical  and 
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sublimation  processes.  For  instance,  sublimation  into  C3,  Eq.  (18),  contributes 

Sea  =  —0-25rjVic3Pw^w,c3j 
ec3  =  Vv,o$  mc3nE,  c3Cw,e3  ■ 

For  a  species  s,  the  total  mass  flow  rate  at  wall  is  the  sum  of  the  mass  flow  rate  in  the  pyrolysis 
gas  and  the  mass  flow  rate  due  to  surface  removal 


Species  flow  rate  at  wall  =  mPt8  +  rh*  a.  (56) 

In  the  gas  phase,  at  an  infinitesimally  small  distance  away  from  the  wall,  the  species  mass  flow  rate 
is  the  sum  of  a  convective  term,  pwp~was tW,  and  a  diffusion  term 


Species  flow  rate  in  gas  phase  =  pwvwa8.w  —  ^  2  )w(-x—)w 

Mp 1  d y 

_  hr/  f  ..  .  ,  /  Ock  s  ,  i 

Pw^W^S  ,W  ^  )w  H”  ^  S,W\—7j^~  )w\ 

Because  expression  (56)  and  expession  (57)  must  be  equal  to  each  other,  one  has 

•  ,  _  Mr  /  ^  ^  /d&s\  -I 

H-  -  Pw^W^SfW  s  ,w )w  +  )w\- 


On  the  other  hand  rht  s  can  be  written  as 


(57) 


(58) 


mt,s 


) 'Yk,s&k,w  H-  rht;,®* 
k 


(59) 


For  example,  jk,s  by  0  +  C(e)  — >■  CO  is 


do,o  —  0.25/^rc)p.„,  fb’,o  ■ 


The  mass  flow  pwvw  can  be  written  in  turn  as 

Pw^w  —  hip  +  ^  ( (  ds  CV3pr  (3)  • 

3 


(60) 


By  substituting  these  into  Eq.  (58),  one  obtains  the  equation  governing  the  wall  value  of  mass 
fractions  of  species  s  as 


[  (rhp  +  'y  )(/h-f>7r..;r;  f  h)f  o  (  r.  )w\&8,w  T  ^  (  t'A'-rr ^A'-rr 

jCy,  oy 

k  k 

P-w  ^  /  dcxs .  _ 

t  n  ) -tt?  —  fflp  s  rnv  s, 

Sc  w  ’  dy  ' 


(61) 


Eq.  (61)  forms  a  system  of  AT  simultaneous  quadratic  equations  in  ckSjU;!s.  To  determine  aw’’ s 
using  this  equation,  one  must  determine  the  wall  slope  ( das/dy)w .  The  wall  slope  is  determined  by 
solving  the  flow  field.  Flow  field  can  be  solved  for  this  purpose  either  by  using  the  so-called  two-layer 
method,  i.e.,  inviscid  layer  plus  boundary  layer,  or  by  using  computational- fluid-dynamics  (CFD). 
In  the  following,  these  two  methods  are  reviewed. 


3.3.  APPLICATION  TO  BOUNDARY  LAYER  ANALYSIS 

In  the  two-layer  method,  the  inviscid  flow  field  is  solved  by  any  available  method.  Then  the 
boundary  layer  is  attached  to  the  inner  boundary  of  the  inviscid  solution.  The  outer  boundary 
condition  of  the  boundary  layer  is  thus  given  by  the  inner  boundary  value  of  the  inviscid  solution.  It 
is  believed  that  several  computer  codes  exist  in  the  world  presently  in  which  this  method  is  used.  To 
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a  high  accuracy,  chemical  equilibrium  can  be  assumed  at  the  edge  of  the  boundary  layer.  However, 
within  the  boundary  layer,  frozen  flow  assumption  is  valid  to  a  quite  high  accuracy. 

Boundary  layer  analysis  was  introduced  in  Section  2.8.  Here,  we  have  to  generalize  it  into 
off-stagnation  points  and  to  the  case  of  Prandtl  number  different  from  1,  and,  above  all,  must  use 
the  bifurcation  model  introduced  above.  Coordinate  transformation  is  unchanged.  The  resulting 
differential  equation  becomes,  for  mass  and  momentum, 


PP  fy 

Pef^e 


+  ff" 


+  2f3(  —  —  fr2)  =  2£(/'^j- 
P 


where 


The  energy  equation  becomes 


d\ln(ue) 

d[ln(Z)} 


/  PP  1 

K  PePe  Pr 


°'v+^,+k  id 


PP'  el  ell  I  / 
PePe  J 


The  species  equation  is  transformed  into 


(62) 


(63) 


PP  Cs  i  ,  (  pp  C8  , 

-«J  +(/+  Tr)aa 


PePe  Sc 


PePe  Sc 


,  PP  C's 
PePe  Sc 


)'  -  2 f%]at 


(64) 


where 

a  _  d[ln{as>e )] 

P'~  d[ln{0]  - 

There  are  N  (number  of  species)  equations  of  the  form  Eq.  (64). 

Eqs.  (62)  to  (64)  form  a  system  of  nonlinear  ordinary  differential  equations  because  a)  Eq. 
(62)  contains  terms  consisting  of  products  of  dependent  variables,  and  (b)  pe/ p  in  Eq.  (62)  depends 
nonlinearly  on  g  and  a’s.  This  system  of  equations  can  be  solved  through  approximate  linearization 
and  iteration.  In  the  end,  at  each  ^-point,  the  solution  of  the  species  equations  can  be  written  as 


as(v)  =  as,wYi(v)  +a'8,w  Y2  (v)  +Yp(v)- 

where  the  homogeneous  solutions  Y% ,  and  the  particular  solution  Yp  are  chosen  so  that 

5f(0)  =  1,  T/(0)  =  0, 


(65) 


y2(0)  =  0,  Yj(0)  =  1, 

yp( o)  =  o,  r;(o)  =  o. 

By  demanding  that  as  value  in  Eq.  (65)  becomes  the  boundary  layer  edge  value  at  the  boundary 
layer  edge,  one  has 

as,toTi(s,776)  +  a'swY2(s,r]e)  +  Yp(s,r]e)  =  aSjC.  (66) 

There  are  N  equations  of  the  form  (66) . 

Eq.  (61)  is  converted  into  the  boundary  layer  coordinates  by  introducing 

das  dr)  , 
dy  dyas 

which  leads  to  a  system  of  simultaneous  nonlinear  algebraic  equations: 


[— (?hy  +  'y^j(SfcOlk,w  +efc)  +  -^—{-^)  w^‘s!w\as,w  +  7 k,sak, 


+  7T- Ck,w{ ^~,)was,w  ~  ~ 


dr) 


Sc. 


dy 


(67) 
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Eqs.  (66)  and  (67)  form  2 N  algebraic  equations  containing  2 N  unknowns  as.,r  and  a'sw.  These 
equations  can  be  solved  through  approximate  linearization  and  iteration.  The  values  of  as  w  are  first 
assumed.  These  values  are  used  in  evaluating  5kCtk,w  Then  there  results  a  system  of  linear  equations 
of  order  2  N 


[A]X  =  R. 


Here, 


For  k  =  1  to  N  :  X  =  (vj. 

For  k  =  N  +  1  to  2 N  :  X  =  a'k  .w. 


This  equation  set  is  solved  as 
to  obtain  improved  values  of  asw 


X  =  [A]^1R, 

and  a's  w.  For  the  next  iteration, 


Next.  as>w  =  e  x  ( New  a8tW)  +  (1  —  e)  x  ( old  a„iW) 


(68) 


where  e  is  a  number  smaller  than  1,  is  used  in  the  evaluation  of  and  so  forth.  This  process  is 

repeated  until  a  converged  set  of  solution  vector  is  obtained.  Experience  has  shown  that  the  most 
appropriate  value  of  e  varies  between  0.2  and  1  depending  on  the  problem. 


3.4.  APPLICATION  TO  CFD 

CED  codes  are  believed  to  exist  worldwide  in  which  the  kinetic  boundary  conditions  cure  applied 
as  outlined  in  this  lecture  (e.g.  Ref.  16).  In  CFD,  the  relation  (61)  becomes  the  boundary  condition 
for  species  variables.  The  wall  slope  ( das/dy)w  is  written  in  a  difference  form.  Because  Eq.  (61)  is 
nonlinear  in  (Xsw,  approximate  linearization  must  be  carried  out.  It  can  be  achieved  by  assuming 
$kOik,w  to  be  known.  In  updating  5kak,w ,  it  would  be  wise  to  correct  only  partly,  i.e.  use  Eq.  (68), 
to  ensure  stability. 

Applying  the  boundary  condition  of  a  linear  form  will  not  be  of  any  difficulty  in  the  CFD 
community.  If  a  two-point  representation  of  wall  slope  is  used,  Eq.  (61)  will  supply  the  diagonal 
term  and  one  off-diagonal  term  in  the  difference  formulation.  If  a  three-point  representation  of  wall 
slope  is  used,  two  off-diagonal  terms  will  be  specified.  If  a  steady-state  solution  is  sought  and  if 
the  solution  is  obtained  by  a  time-asymptotic  method,  the  approximate  linearization  scheme  will  be 
implemented  by  updating  SkakjW  after  each  iteration.  If  a  time-accurate  solution  is  sought,  a  unique 
iterative  procedure  must  be  devised  in  which  the  iteration  is  performed  while  time  is  frozen. 


4.0.  SURFACE  PROBLEMS  IN  PLANETS 
4.1.  EARTH  ENTRIES 

In  Earth’s  atmosphere,  the  gas-surface  interaction  problem  was  studied  in  the  past  mostly  in 
two  areas:  catalytic  recombination  of  oxygen  and  ablation  of  carbonaceous  heatshield.  Goulard’s 
pioneering  work  for  the  stagnation  point  was  discussed  already.  His  work  was  extended  to  off- 
stagnation  points,  three  dimensional  flows,  and  to  nitrogen  recombination  by  various  researchers. 

Until  recently,  ablation  phenomenon  was  treated  by  invoking  gas-surface  equilibrium.  The  as¬ 
sumption  of  gas-surface  equilibrium  circumvents  the  need  to  know  the  //,.  /  and  In  the  kinetic 

approach,  the  equilibrium  condition  can  be  reached  by  assuming  r)rj  and  rjvj  to  be  a  large  number. 
According  to  Eqs.  (8a),  //,../  and  rjvj  need  only  to  be  much  larger  than  Mm / \/ Rem  to  produce  an 
effect  of  a  gas-surface  equilibrium.  The  equilibrium  approach  and  the  kinetic  approach  have  been 
compared  in  Ref.  13.  It  found  that  equilibrium  approach  gives  always  higher  ablation  rate  than  the 
kinetic  approach. 

The  recent  data  on  nitridation10 ,  N  +  C(c)  — >■  CM,  which  gives  rjrj  to  be  0.3,  has  brought  an 
alarm  and  skepticism  of  the  measurement13.  This  is  because,  if  the  measured  value  of  rjrj  is  accurate, 
this  process  will  be  the  major  process  of  ablation  in  a  hyperbolic  Earth  entry.  More  experiment  is 
desirable. 

For  wall  heat  transfer  rates  in  the  range  from  5  x10s  and  1.5xl06  W  /nr.  such  as  in  the  stagnation 
region  of  the  nose  cone  of  the  Space  Shuttle,  carbon-carbon  composite  overcoated  with  silicon  carbide, 
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SiC,  is  used  for  the  surface.  When  exposed  to  a  stream  containing  atomic  oxygen,  SiC  is  converted 
quickly  to  Si02-  Si02  is  in  an  amorphous  form.  Oxygen  atoms  are  dissolved  in  this  amorphous  Si02, 
diffuse,  and  reach  the  SiC  substrate.  Once  they  reach  the  SiC  substrate,  they  combine  with  C  to 
form  gaseous  CO.  These  CO  molecules  are  dissolved  in  Si02,  diffuse,  and  escape  from  the  surface, 
leaving  Si02-  This  phenomenon  has  been  studied  experimentally  to  a  considerable  extent  in  Europe 
and  in  Japan.  Yet,  there  are  no  theoretical  model  to  describe  the  rate  process  quantitatively. 

In  the  future,  attempts  will  be  made  to  construct  an  airbreathing  hypersonic  vehicle.  High 
heating  rates  will  occur  at  the  nose  tip,  wing  leading  edge,  cowl  lip,  and  in  the  combustion  chamber 
of  such  a  vehicle.  New  material  such  as  halfnium  di-boride,  HfB2,  may  be  used  for  such  high  heating 
areas.  Gas-surface  interaction  for  such  a  new  material  will  be  a  challenge. 

The  ablation  rate  of  carbonaceous  heat. shield  material  in  the  nose  region  of  a  blunt  body  was 
measured  systematically  in  Ref.  17.  The  measured  ablation  rate  was  about  3  times  higher  than 
calculated  assuming  gas-surface  equilibrium  in  the  stagnation  region.  This  cannot  be  attributed  to 
roughness-caused  turbulence  because  surface  roughness  cannot  produce  turbulence  at  the  stagnation 
point,  because  velocity  is  zero  there.  In  Ref.  18,  this  large  ablation  rate  is  attributed  to  the  turbulence 
generated  by  ablation. 

Various  experiments  showed  that  the  highest  heating  rates  occur  at  points  away  from  the  stag¬ 
nation  point,  somewhere  where  the  surface  normal  is  about  40°  with  respect  to  the  flow  direction 
(see,  e.g,  Ref.  19).  This  peak- heating  region  is  usually  the  point  where  turbulent  transition  occurs 
due  to  roughness.  In  many  cases  it  coincides  also  with  the  point  of  sphere-cone  juncture  where  an 
adverse  pressure  gradient  occurs.  Heat  transfer  rate  in  the  region  downstream  of  this  point  is  usually 
higher  than  the  heat  transfer  rate  to  the  stagnation  point,  and  therefore  the  ablation  rate  is  also 
higher  there.  Theoretical  prediction  of  the  ablation  rate  in  this  downstream  region  from  the  first 
principles  is  at  this  time  not  yet  achieved.  To  do  so,  one  needs  to  correctly  model  injection- induced 
turbulence  and  roughness-induced  turbulence. 

If  one  can  model  turbulence  accounting  for  these  effects,  then  the  theory  given  in  Chapter  3  can 
be  applied.  The  only  part  that  needs  modification  is  viscosity  p  which  will  have  to  be  modified  to 
include  the  effect  of  turbulence. 

Because  there  is  a  large  uncertainty  about  the  behavior  of  turbulence  in  the  downstream  region, 
the  Ik shield  is  usually  over-designed.  In  the  past,  this  philosophy  was  proven  to  be  effective  and 
necessary. 

4.2.  MARTIAN  ENTRIES 

Martian  entries  occur  typically  at  6  krri/s.  The  density  of  the  atmosphere  is  about  1/70  of 
that  of  the  Earth.  Heating  rate  does  not  exceed  4xl06  W/m2.  A  material  made  of  silica  matrix 
impregnated  with  an  organic  resin  is  used  as  an  ablating  heatshield.  As  mentioned  earlier,  Martian 
atmosphere  consists  mostly  of  CO2.  At  the  edge  of  boundary  layer,  CO2  is  dissociated  into  CO  and 
O.  CO  and  O  diffuse  through  the  boundary  layer  and  reach  the  wall.  At  wall,  CO  +  O  —»■  CO2  can 
conceivably  occur.  The  surface  of  the  material  during  entry  is  pure  silica.  As  mentioned,  silica  shows 
a  near- zero  reaction  coefficient  for  CO  +  O. 

One  uncertainty  is  about  the  effect  of  surface  roughness  on  the  effective  recombination  coefficient, 
Eq.  (21).  Because  we  do  not  know  the  BET  area  of  a  Martian  heatshield  during  an  entry  flight, 
we  cannot  determine  the  effective  recombination  coefficient.  In  addition,  we  do  not  know  where 
turbulent  transition  occurs  because  we  do  not  know  the  intensity  of  the  injection-induced  turbulence 
or  the  roughness-induced  turbulence. 

4.3.  VENUSIAN  ENTRIES 

The  atmospheric  composition  of  the  planet  Venus  is  similar  to  that  of  Mars.  But  the  density 
of  the  atmosphere  is  about  70  times  that  of  the  Earth,  and  the  entry  speed  is  more  than  10  km/s. 
At  the  peak-heating  point,  the  pressure  at  the  stagnation  point  is  of  the  order  of  10  atm.  The 
combination  of  the  density  and  entry  speed  produces  heating  rates  in  the  range  of  several  million 
W .in2.  Only  carbonaceous  ablation  material  can  be  used  as  heatshield.  The  comments  about  carbon 
surface  made  in  Chapters  1  and  2  are  applicable. 

In  the  Pioneer- Venus  Probe  mission,  four  probe  vehicles  entered  into  the  Venusian  atmosphere. 
Each  vehicles  were  equipped  with  two  thermocouples  in  their  heatshield,  one  near  the  stagnation 
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point  and  the  other  near  the  frustum  edge.  The  data  from  these  thermocouples  indicated  that  the 
heat  transfer  rates  were  higher  at  the  frustum  edge  than  at  the  stagnation  point.  The  flight  data 
are  consistent  with  the  assumption  of  a  turbulent  flow.20,21 

In  this  entry  flights,  the  edge  of  boundary  layer  will  contain  oxygen  atoms  and  a  small  concen¬ 
tration  of  nitrogen  atoms.  The  reactions  (13)  and  (14)  will  occur.  Reaction  (13)  will  contribute 
considerably  to  the  total  heat  transfer  rate.  Sublimation  reaction  (18)  will  occur  also. 

4.4.  TITAN  ENTRIES 

The  atmosphere  of  Titan  consists  mostly  of  N2,  with  a  small  amount  of  GH4.  The  density  of  the 
atmosphere  is  nearly  the  same  as  that  of  the  Earth.  For  the  Huygens  vehicle  sent  by  the  European 
Space  Agency,  the  entry  velocity  was  about  5  km/s.  However,  in  future  missions,  entry  velocity  of 
up  to  about  14  km  is  possible.  Radiation  from  the  CN  becomes  a  major  component  of  heat  transfer 
rate.  Both  a  silica-based  and  a  carbon-based  heatsliield  can  be  used  for  these  entries.  There  will 
be  a  small  amount  of  dissociation  of  Ng,  and  full  dissociation  of  CH4.  Reactions  (15)  and  (16)  will 
occur  in  both  directions. 

4.5.  OUTER  PLANET  ENTRIES 

Atmosphere  of  the  outer  planets  consists  of  hydrogen  and  helium.  Entry  into  Jupiter  has  been 
made  in  the  Galileo  Probe  mission.  Heat  transfer  was  due  overwhelmingly  by  radiation.  The 
heatsliield  was  made  with  a  carbon-based  material.  During  this  entry  flight,  progression  of  surface 
recession  was  measured.  The  flight  data  showed  a  surprisingly  low  surface  recession  at  the  stagnation 
region  and  surprisingly  high  recession  at  the  frustum. 

The  low  recession  in  the  stagnatin  region  is  attributed  to  thermal  nonequilibrium.22  The  high 
recession  in  the  frustum  region  is  attributed  to  increase  in  radiation  due  to  turbulence  caused  by 
ablation.23  The  dominant  gas-surface  interaction  was  sublimation,  reaction  (18).  The  sublimation 
rate  was  driven  essentially  by  the  energy  conservation  law,  rather  than  by  surface  kinetics. 

In  a  futuristic  mission,  one  can  conceive  entry  flights  into  Saturn,  Uranus,  or  Neptune.  In  such 
entries,  radiative  heating  will  be  much  less  than  that  to  the  Galileo  Probe.  There,  not  only  the 
sublimation  but  the  eactions,  Eqs.  (15)  and  (16)  could  conceivably  become  important.  Surface 
kinetics  may  play  a  role. 

5.0.  CONCLUDING  REMARKS 

In  the  above,  gas-surface  interaction  problem  was  studied  from  the  kinetic  approach.  The  kinetic 
approach  enables  rational  way  of  quantitatively  evaluating  the  surface  phenomena  in  determining 
heat  transfer  rates  and  surface  ablation  rates.  The  study  identifies  areas  of  uncertainty  also,  which 
are:  1)  the  effective  coefficients  of  gas-surface  interactions  for  rough  surfaces,  and  2)  the  intensity  of 
turbulence  for  a  rough  surface  in  the  presence  of  ablation.  As  a  minor  and  overcomable  difficulty, 
nonlinearity  in  the  governing  equations  at  wall  is  pointed  out. 

The  problem  of  gas-surface  interaction  in  hypersonic  flight  contains  several  different  disciplines 
that  are  tied  together  in  a  rather  unexpected  way.  Because  of  the  unusual  nature  of  the  problem,  the 
traditional  discipline-based  approach  to  the  problem  becomes  difficult  to  apply.  This  is  probably  the 
reason  for  the  misunderstanding  and  mismanagement  of  the  issues  that  existed  in  the  past.  To  ensure 
mission  success  in  the  midst  of  such  uncertainties,  the  planetary  entry  vehicles  had  been  designed  in 
the  past  with  a  rather  large  safety  factor.  In  retrospect,  we  feel  that  such  large  safety  factors  were 
justified,  because  several  unexpected  phenomena  occurred.  For  a  rational  management  of  the  future 
space  programs,  all  scientific  issues  about  gas-surface  interaction  will  have  to  be  understood  better 
by  all  involved.  Much  more  research  on  the  subject  is  highly  desirable. 
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